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Abstract 

It is assumed that axions exist and are created in the sun by Primakoff conversion 
of photons in the Coulomb fields of nuclei. Detection rates are calculated in germanium 
detectors due to the coherent conversion of axions to photons in the lattice when the 
incident angle fulfills the Bragg condition for a given crystalline plane. The rates are 
correlated with the relative positions of the sun and detector yielding a characteristic 
recognizable sub-diurnal temporal pattern. A major experiment is proposed based on a 
large detector array. 

Introduction 

The putative pseudo-goldstone axion [1,2] had its inception almost twenty years 
ago following the suggested solution to the strong CP problem by Peccei and Quinn [3]. 
This prompted many theoretical investigations and experimental searches. 

An extensive review of axion phenomenology, and their effects on stellar 
processes, and bounds, was given by Raffelt [4]. A detailed treatment of solar axions, and 
of a proposed method of detecting them, was given by van Bibber et al. [5]. In this letter 
we analyze a completely different direct-detection technique based on coherent Primakoff 
conversion of axions to photons in the lattice of a germanium detector when the incident 



angle satisfies the Bragg condition. Some of the basic concepts have been presented 



earlier in different contexts [6,7]. 

The axion-photon interaction Lagrangian is Li n t = (1/4 M) a F^ v F a P G^ V a(3 ' where 

"a" is the pseudoscalar axion field, F^ v is the electromagnetic field, and 1/M = g a yy is the 
axion-photon coupling. In the Born approximation the differential cross section for 
elastic axion-to-photon conversion off an atom of charge Z (see Fig. 1) is: 
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where 20 is the scattering angle and F(20) is the atomic form factor. We have 
independently verified this cross section given previously [6,7]. Note that the axion-to- 
photon cross section is twice as large as that for photon-to-axion conversion because the 
photon has two helicities. 

The form factor, F a (20), is that used in earlier work [7] and can be written in two 

forms: 



F a (20) = 
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= F a (q,k), 
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where ro is the screening length, and q=2k sin is the momentum exchange. The 
differential cross section can then be rewritten as follows: 
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Averaging equation (3) over all solid angle yields: 
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where r| = r k, and Z 2 ah 2 /SnM 2 c 2 = 1.15 x 10" 44 cm 2 = VTg for Z = 32, where 
V^= (gayyX 10 8 GeV) 2 is the (unknown) relative strength of the axion-photon coupling. 



The cross section for X = 1 integrated over the energy range to 12 keV is 
1.46 x 10- 43 cm 2 . 

The expected flux of axions from the sun has been calculated by van Bibber et al. 
[5]. The spectrum is continuous with a peak near 4 keV and falls off exponentially 
beyond about 10 keV. The results of these detailed numerical calculations are well 
approximated by the empirical form: 



d4> = /^$o ( E/E o) 
dE 



E e E/Eo _ j 



(5) 



where 3>o = 5.95 x 10 14 cm -2 sec -1 . We have used this same expression but have adjusted 
the parameter Eo = 1.103 keV to take into account a small change in the solar model core 
temperature when helium and metal diffusion are included [8]. 

For very light solar axions the Primakoff process in a periodic lattice is coherent, 
similar to Bragg reflection of x-rays, which leads to the Bragg condition that the 

momentum transferred to the crystal must be a reciprocal lattice vector, G = — f h, k, l) 

a 

where ao is the size of the conventional cubic cell and h, k, and 1 are integers [9]. If we 
assume the charge distribution in the crystal is simply a superposition of atomic charge 

distributions, the rate of conversion per unit flux of axions of energy E a and momentum 

— > -> 
k a to photons of energy Ey and momentum k y is: 
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where V is the volume of the crystal, v c is the volume of the unit cell, and S(G) is the 
structure function for the crystal, which for germanium is: 
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Table 1 gives the parameters of the most relevant reciprocal lattice vectors of germanium. 

If we integrate over all final photon states we find the total conversion rate of 
axions with momenta k^ and energy E a to be 
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Note that in equation (8) the conversion rate is essentially zero unless the energy and 
momentum of the axion exactly satisfies the Bragg condition. In a real crystal, the 
outgoing x-ray very probably produces a photo-electron after traveling a distance of order 
microns which tends to smear out the delta function. However, so long as the mean free 
path of the x-ray is long compared with the size of the unit cell (as it is in germanium), the 
broadening of the delta function will always be much less that the actual energy resolution 
of the detector and can be neglected. 

If we now multiply this by the flux of solar axions, (5), the total rate of 
conversion of axions is, 
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where k is the unit vector from the sun, and d<£/dE is evaluated at the axion energy 



-> 



E a = he |G| 2 / 2k • G . This can be written in the compact form: 
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where is the mass of the detector in kg, g = (h, k, 1), y = ao/ro, p = — ^— 

agEo 
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, and Nq = — - — ^-y = 0.61 / kg ■ d for germanium. 
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For a given direction of incident axions, k, the Bragg condition for each reciprocal 

A -> 

lattice vector such that k • G > determines a narrow range of energies. A germanium 
detector has a characteristic low-energy resolution, AE ~ 400 eV. We take this into 
account by replacing the delta function in (10) by a gaussian with the same fwhm as the 
detector. Finally we calculate the total counting rate in a range of energies of width Ae , 



e+Ae/2 

R(k,e)= jd e ' dE (k, G '). (11) 
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In figure 2 we show the expected counting rate, R (k(t),e ) = R(t), as a function of 
time over a single day. The position of the sun is calculated using the U. S. Naval 
Observatory Vector Astronomy Subroutines (NOVAS) [10]. 

The pronounced variation in the axion counting rate as a function of time suggests 
that events in the detector be analyzed in terms of the correlation function: 

% = | o T [R(t)-R]n(t)dt, (12) 

where R(t) is the theoretical instantaneous axion counting rate in a given energy interval, T 
is the time during which the detector is on, Ris the average of R(t) over this period and 
n(t) is the number of counts at time t in a short time interval At. Note that if n(t) is 
uncorrelated with the position of the sun, then % ~ within statistical fluctuations, 
whereas if n(t) contains an axion component, % will increase proportionally to T. 

We have carried out Monte Carlo simulations of a germanium detector with 
realistic background rates taking the range of energies from 2 to 10 keV in 0.5 keV 
intervals. We find that an array of 400, oriented single crystal n-type detectors, each of 
0.25 kg operating for 10 years can set an upper limit on the scale factor X, 

X<10" 4 . (13) 
This corresponds to a 95% CL bound on the axion-photon coupling constant 
g a yy = 1/Mc 2 ~ 5 x 10" 10 GeV" 1 for axion masses between and about 1 keV. 

Similar results could be obtained using -160, single crystals of p-type detectors, 
each of 0.65 kg. A negative result would exclude a significant portion of the remaining 
g a yy - m a phase space. A pilot experiment has been done which has accumulated data for 
approximately 1.92 kg-years and is being analyzed by the methods discussed above [11]. 
This result was derived using the correlation function of equation (1 1). It is not clear that 
more sensitive analysis techniques employing pattern recognition, for example, might not 
be more effective, these are under investigation at present. 

This large array of low energy threshold, low background detectors with excellent 
energy resolution would also be a sensitive detector for weakly interacting massive 
particles hypothesized as cosmic cold dark matter. 



The present technique could also be applied to low temperature sapphire 
detectors planned for dark matter searches, and to tellurium oxide detector arrays being 
constructed for double-beta and dark matter searches. In these searches, nuclear recoil 
events are of interest, whereas for an axion search, only events associated with photon 
signals are relevant. The same is true for large arrays of Ge detectors which also employ 
cryogenic bolemetric techniques. For these to be more effective than an ordinary low- 
background Ge detectors, the raw background due to photons would have to be low. 
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Captions for the Figures 



Figure 1. Axion conversion to a photon in the Coulomb field of the nucleus by the 
Primakoff effect. 

Figure 2. Theoretical prediction of the count rate of photons converted from axions 
incident at a Bragg angle, for a detector located at Sierra Grande, Argentina 
(41« 4i' S , 65° 22' W). The rate was calculated for 1/Mc 2 = g a yy = 10" 8 
GeV -1 . The location that was chosen is where the pilot experiment is being 
performed. 
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